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Abstract 

Using Renyi entropy, a possible thermostatistics for nonextensive 
systems is discussed. We show that it is possible to have the q- 
exponential distribution function for nonextensive systems having non- 
additive energy but additive entropy. It is also shown that additive 
energy, as a approximation within nonextensive statistics, is not suit- 
able for discussing fundamental problems for interacting systems. 
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1 Introduction 

Renyi entropy [1] 

In V" r)'^ 

^R^^_nU=iP^^ g>0 (1) 
1 - q 

often applied in the study of multi-fractal and chaotic systems[2], where Pi 
is the probability that the system is at the state labelled by i (Boltzmann 
constant k=l) and w is the total number of the states. In these studies, 
Renyi entropy was associated with the exponential probability distributions 
of BGS[2] which as a matter of fact are not the distributions derived from 
this entropy. Since the proposal [3] of the nonextensive statistical mechanics 
(NSM), there is a growing interest in Renyi entropy which has been com- 

pared to Tsallis entropy [3] 5*"^ = — , associated with a g-exponential 

distribution expq{x) = [1 + (1 - g)a;]i/(^-«), in the discussions of possible 
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nonextensive statistics and the relative fundamental problems such as ther- 
modynamic equilibrium and stability[4, 5, 6, 7, 8, 9, 10, 11, 12] for systems 
having additive energy or infinite number of states[6]. 

Regarding the possible statistics of Renyi entropy, some questions should 
be asked : is additive just as, e.g., the entropies of Boltzmann-Gibbs- 
Shannon statistics (BGS) for systems having product joint probability [13], 
but should it be associated to independent systems having additive energy 
just as in BGS? Should it be associated to exponential probability distribu- 
tions as has been done in many works? Are there alternative distributions 
intrinsic to Renyi entropy if it can be maximized to get thermodynamic 
equilibrium? If the Renyi statistics is not extensive, what are its possible 
nonextensive properties? Due to the importance of in the study of chaos 
and fractals and since is identical to Boltzmann entropy S = \nW for the 
fundamental microcanonical ensemble[2], the possible answers to the above 
questions would be interesting for both BGS and its possible extensions based 
on S^. 

In this paper, I will present a thermo-statistics derived from Renyi entropy 
for equilibrium systems with nonadditive energy required by the existence of 
thermodynamic equilibrium[l i]. This formalism may be considered as an 
alternative to BGS and to NSM for canonical systems with additive entropy 
but nonadditive energy [15]. 



2 Renyi and Tsallis entropies 

There is a monotonic relationship between these two entropies : 

1 — q ^ — Q 

and, for complete probability distribution {J2T=iPi — 1) microcanonical 
ensemble, is identical to the Boltzmann entropy S : 

S^ = S = \nw (3) 

since J2T=iPi = w^~'^- Other properties of can be found in [1, 2, 4]. 

The concavity of and S'^ for g > 1 is shown in Figure 1 and 2, 
respectively (they are convex and have minimums for g < 0, so we do not 
consider this case in this paper). It is worth noticing that the two entropies 
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get their maximum at the same time for any q. But the maximum of S is 
In W and independent of g, while the maximum of S'^ is and decreases 

down to zero when g — oo and increases up to — 1 when g — > 0. It should 
be mentioned that is not always concave for g > 1, as shown in Figure 1. 
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Figure 1: The concavity of Renyi entropy S for g > with a two probability 
distribution pi = x and p2 = ^ — x. Note that the maximal value does not 
change with g. The maximum becomes more sharp for larger g. The curve 
for g > 1 shows that is not always concave, but the maximum remains 
the unique extremum. 

Due to the fact that is a monotonically increasing function of S""^, 
they reach the extremum together (see Figures 1 and 2). One can hope 
that the maximum entropy (for g > 0[4]) will give same results with same 
constraints. Indeed, Renyi entropy has been used to derive, by maximum 
entropy method, the Tsallis g-exponential distribution within an additive 
energy formalism[16, 17]. In what follows, we present a thermostatistics 
based on for nonextensive systems having additive entropy and nonaddi- 
tive energy.. 



3 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 



Figure 2: The concavity of Tsallis entropy S for g > with a two probabihty 
distribution pi = x and p2 = 1 — x. Note that the maximal value increases 
from zero to unity with q decreasing from infinity to zero. 

3 Canonical distribution of Renyi 

We suppose complete distribution J2T=iPi = 1 ^.nd U = J2T=iPi^i where U 
is the internal energy and Ei the energy of the system at the state i. We will 
maximize as usual the following functional : 

F='^^^ + aY.P^'lT.P^E,- (4) 

? i=i i=i 

We get 

PiO,[a-^E:]^l^''-^\ (5) 

Since recovers Boltzmann-Gibbs entropy S = — Z^lLiPilnpi when q = 1, 
it is logical for us to require that Eq.(5) recovers the conventional exponential 
distribution for g = 1. This leads to 

p^ = ^[l-iq-l)PE,Y/^'^-'\ (6) 
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where {q-l)p = 7/a and Z = ET=i[^ -(q- l)pE,Y/^''-^\ We will show the 
physical meaning of (3 later. 

r, 2 q 

Note that the second derivative ^-5- = — r ^' g [1 + -r— r ^' ■? ] is negative for 
any distribution only when g < 1. In general may be positive for g > 1 

so that the above distribution Eq.(6) is not stable. 

This formalism can be applied to an important family of chaotic behaviors 
described by Levy flight[l(S] if we assume < g < 2/3. In this case, the first 
moment x calculated with the usual expectation is convergent for 1/2 < 
g < 2/3 and can be used as a constraint of maximum entropy to obtain 
p{x) ~ [1 + (1 — q)(3x]~^^^^^^\ For Levy flight distribution with large x, 
p ~ x"'^^'^, one obtain : g = 7/(1 + 7) (0 < 7 < 2). 

4 Mixte character : nonadditive energy and 
additive entropy 

It has been shown[l 1] that, for thermal equilibrium to take place in nonexten- 
sive systems, the internal energy of the composite system A + B containing 
two subsystems A and B must satisfy 

U{A + B) = U{A) + U{B) + XU{A)U{B) (7) 

which means 

Eij{A + B)= E,{A) + Ej{B) + \Ei{A)Ej{B) (8) 

where A is a constant. Applying Eq.(8) to Eq.(6), we straightforwardly get 
the product joint probability : 

Pi,{A + B)=p,iA)p^{B) (9) 

and the additivity of : 

S'^iA + B) = + (10) 

if A = (1 — q)j3. So is essentially different from 5"^, because in this case 
is nonadditive with S^{A + B) = S^iA) + S^{B) + (1 - q)S^{A)S^{B). 
Note that we do not need independent or noninteracting or weakly interacting 
subsystems for establishing the additivity of or the nonadditivity of 5*"^, as 
discussed in [14, 19, 20]). So in this formalism, we can deal with interacting 
systems with nonadditive energy but additive entropy. 
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5 Zeroth law and temperature 



It is easy to show that, from Eq.(6), 

^^ = lnZ + ln[l + (l-g)/5f/]/(l-g). (11) 



So we have 

P=[l + {l-q)m^ (12) 

or 

i^^-d-,)^'. (13) 

Since [1 + (1 — q)PU] is always positive (g-exponential probabihty cutoff), 
has always the same sign as j3 can be proved to be the effective inverse 
temperature if we consider the zeroth law of thermodynamics. Let 5S^ be 
a small change of of the isolated composite system A + B. Equilibrium 
means 5S^ = 0. From Eq.(lO), we have SS^{A) = —6S^{B). However, 
from Eq.(7), the energy conservation of A + B gives SU{A) = 

- [iHi-^mm ^^^^^- That leads to 

[1 + (1 - qWiA)]^^ = [1 + (1 - ^)mB)f-l^ (14) 
or j3{A) = (5{B) which characterizes the thermal equilibrium. 



6 Some "additive" thermodynamic relations 

Due to the mixte character of this formalism with additive entropy and non- 
additive energy, all the thermodynamic relations become nonlinear. In what 
follows, we will try to simplify this formal system and to give a linear form 
to this nonlinearity. 

Using the same machinery as in [20] which gives an extensive form to the 
nonextensive Tsallis statistics, we define an additive deformed energy E as 
follows : 

E = \n[l + {l-q)l3U]/{l-q)l3 (15) 

which is identical to U whenever q = l. Note that E{A + B) = E{A) + E{B). 
In this way, Eq.(ll) can be recast into 

= In Z + pE. (16) 
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So that [3 = The first law can be written as 

6E = T6S^ + Y6X (17) 

where Y is the deformed pressure and X the coordinates (volume, surface 
...) and T = The deformed free energy can be defined by 

F = E- TS^ = -T In Z, (18) 

so that Y = The real pressure is Y^ = [1 + {1 — q)(3U]Y and the work 
is 6W = Y^6X. The deformed heat is 6Q = T6S^ and the real heat is 
(5g« = [1 + (1 - q)(3U]6Q 

7 Grand-canonical distributions 

It is easy to get the grand-canonical ensemble distribution given by 

P^ = ^[l-{q- - ^iN,)Y'^'~'\ (19) 

which gives 

S^ = \nZ + ln[l + (1 - q)l3U]/{l - q) + ln[l - (1 - q)puN]/{l - q). (20) 

Let M be the deformed particle number : M = ln[l — (1 — g)/3u;iV]/ {l — q)j3oj, 
Eq.(20) becomes 

= \nZ + pE + piuM. (21) 
M must be additive, so that N is nonadditive satisfying 

N{A + B) = N{A) + N{B) + (1 - q)/3ujN{A)N{B) (22) 

Due to the distribution function of Eq.(19), the quantum distributions will 
be identical to those in NSM[21]. 

8 Renyi statistics with additive energy? 

Now we know that should be intrinsically associated with the g-exponential 
distributions. Only when g = 1 this statistics recovers BGS and the q- 
exponential becomes the usual exponential function. Then an interesting 
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question is whether or not this statistics may be associated with additive 
energy when q^\. 

Supposing A and B are independent, i.e., Eij{A + B) = Ei{A) + Ej{B), 
let us see the probabihty of the system A + i? for a joint state ij : 

P^AA + B) = ^^^^ [1 - (g - mEU) + E,{B))YI^^~'^ (23) 
= pi{A)pj\i{B I A) 

where 

P^{A) = ^[1 - (g - l)(3EU)Y'^'~'^ (24) 

is the probabihty for A to be at the state i and 

P,\^{B I A) = ^^^^^^^^ [1 - (g - l)/3e,|,(A | S)]^/^'^"^) (25) 

is a conditional probability for B to be at a state j with energy ej\i{A \ B) = 
Ej{B)/[l - {q- l)(3Ei{A)] if yl is at i with energy Ei{A). In this case, the 
total entropy is given by 

S-(A + B) = '"E.P.WE.Kl.t^l-^)'! (26) 

^ S^iA) + S^iB). 

This is in contradiction with Eq.(lO). So is no more additive with additive 
energy. As a matter of fact, Eq.(26) is wrong because Eq.(23) does not hold 
if we consider the product probability Eq.(9). We would get 

P^B) = -^[l-(g- !)/?£;, (S)]V(.-i) (27) 

= z,iB\A) ^^ " ^'^ " ^^^'^^'^"^ ' ^^^^ 

which implies Ej{B) = e^^^A \ B) = Ej{B)/[l - (g - l)(3Ei{A)\ and is 
valid only when g = 1. This means that additive energy may force back the 
nonextensive statistics to BGS if the product joint probability applies. 

We would like to discuss in passing the problem of thermodynamic insta- 
bility of Renyi entropy which has been shown [(i] to be non-observable and 
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instable because an arbitrarily small variation S in probability distribution 
may lead to an important variation in for a system having infinite num- 
ber of states w. According to this analysis, Renyi entropy can be physically 
useful for finite systems having finite number of states. 

However, it should be noted that the Lesche's conclusion for infinite w is 
reached with the asymptotic behavior of for the case oi w ^ oo and of 
finite variation 6 of probability distributions, i.e., 1/w is negligible compared 
to 6. This is a very harsh condition if we consider that S must be arbitrarily 
small for observability condition. It should be noted that the asymptotic 
behavior of for finite S and w — ^> oo is different from the asymptote 
for arbitrarily small S and arbitrary w which can be very large, e.g., 6 is 
smaller than or of same order of magnitude as 1/w. This second asymptotic 
behavior should be more general to our opinion because it applies for any 
system. Taking the probability distributions proposed by Lesche and making 
the same calculations without neglecting anything, one gets, for both g > 1 
and g < 1, AS^{6,w)/ Smax oc {5/2Y for arbitrarily small 5^0. The 
observability condition [(i] is satisfied. This result is in addition consistent 
with the fact that 5*^ is a monotonic function of S'^ which is observable 
according to the same analysis [7]. We indeed have dS^ = ^^(f^^^^T = 

So if dS'^/S'^ 0, we also have dS^ / for finite and S'^. 



9 Conclusion 

In summary, the additive Renyi entropy is associated with nonadditive en- 
ergy to give an nonextensive thermostatistics characterized by g-exponential 
distributions which have been proved to be useful for many systems show- 
ing non Gaussian and power law distributions [22]. This formalism would be 
useful for interacting nonextensive systems whose information and entropy 
may be additive or approximately additive. The problem of the instability 
of Renyi entropy is reviewed. We think that this entropy may be physically 
useful for any number of states. 

A important point should be underlined following the result of the present 
work. Renyi entropy is identical to Boltzmann one for microcanonical ensem- 
ble. So, theoretically, its applicability to systems with nonadditive energy 
means that Boltzmann entropy may also be applied to nonextensive micro- 
canonical systems as indicated by Gross[15]. 
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